8G-FM-39 


JSC-16758 


AUG 7 1980 


Orbit 

Error Analysis 


(NAS A-TM-31030) ORBIT INU ALIGNMENT: ERROR N80-23415 

ANALYSIS (NASA) 42 p HC A03/MF A0 1 CSCL 20A 


Uaclas 

G3/17 26922 


} 


Mission Planning and Analysis Division 
August 1980 


NASA 

National Aeronautics and 
Space Administration 

Lyndon B. Johnson Space Center 

Houston, Texas 





n 


80-FM-39 JSC- i 6758 

SHUTTLE PROGRAM 

t]k 1 ft // 1 

ORBIT IMU AfcWattNI: ERROR ANALYSIS 

By R. W. Corson, McDonnell Douglas Technical Services Co. 

JSC Task Monitor: T. J. Blucker, Mathematical Physics Branch 


Approved:. 


Emil Schiesser, Chief 
Mathematical Physics Branch 


Approved : 



Mission Planning and Analysis Division 
National Aeronautics and Space Administration 
Lyndon B. Johnson Space Center 
Houston, Texas 
August 1980 



CONTENTS 


Section Page 

1.0 SUHMABI 1 

2.0 INTRODUCTION 2 

3.Q DISCUSSION 4 

4.0 RESULTS 6 

4.1 THE RMS IMU ALIGNMENT . ERROR 6 

4.2 THE SINGLE AXIS IMU ALIGNMENT ERROR 8 

4.3 THE MAXIMUM EXPECTED IMU ALIGNMENT ERROR 12 

4.4 THE STS-1 CASE 13 

4.5 SPECIAL CASES^ 13 

5.0 CONCLUSIONS AND RECOMMENDATIONS 16 

5.1 CREW PROCEDURES 16 

6.0 REFERENCES 18 

APPENDIX A - STAR TRACKER, IMU, AND NAV BASE ERROR MODEL PARAMETERS , . 19 

APPENDIX B - MONTE CARLO SIMULATION DERIVED IMU ALIGNMENT ERROR 

STATISTICS 21 

APPENDIX C - TOTAL AND SINGLE-AXIS IMU ALIGNMENT ERROR PROBABILITY 

DENSITY FUNCTIONS FOR STS-1 TYPE ALIGNMENTS 28 

APPENDIX D - DERIVATION OF THE RMS IMU ALIGNMENT ERROR INDICATOR ... 33 

D.O THE RMS IMU, ALIGNMENT ERROR 34 

D. 1 DERIVATION OF ALIGNMENT ERROR INDICATOR 34 

D.2 THE STS-1 CASE 36 


ill 


PRECEDING PAG £ Cu 


l 'w* I I 



TABLES 


Table Page 

4.0- 1 IMU ALIGNMENT ACCURACY INDICATORS 6 

4.1- 1 RMS TOTAL IMU ALIGNMENT ERROR 7 

4.5-1 RMS ALIGNMENT ERROR - SPECIAL CASES 15 

B-1 KEY TO SIMULATION PROGRAM OUTPUTS 22 

B-2 IMU ALIGNMENT ERROR STATISTICS FOR 90° STAR PAIR ANGULAR 

SEPARATIONS 23 

B-3 IMU ALIGNMENT ERROR STATISTICS FOR 75'' STAR PAIR ANGULAR 

SEPARATIONS . . ■ 24 

8-4 IMU ALIGNMENT ERROR STATISTICS FOR 60° STAR PAIR ANGULAR 

SEPARATIONS 25 

B-5 IMU ALIGNMENT ERROR STATISTICS FOR 45° STAR PAIR ANGULAR 

SEPARATIONS 26 

B-6 IMU ALIGNMENT ERROR STATISTICS FOR 30° STAR PAIR ANGULAR 

SEPARATIONS 27 


lv 



FIGURES 


Figure Page 

2.0- 1 IMU ALIGNMENT COORDINATE SYSTEMS AND TRANSFORMATIONS .... 3 

4.1- 1 RMS IMU ALIGNMENT ERROR 9 

4.2- 1 STANDARD DEVIATION OF THE SINGLE AXIS IMU ALIGNMENT ERROR . . 11 

4.3- 1 ESTIMATED MAXIMUM TOTAL IMU ALIGNMENT ERROR 14 

C-1 TOTAL IMU ALIGNMENT ERROR PROBABILITY DENSITY FUNCITCN ... 29 

C-2 X-AXIS IMU ALIGNMENT ERROR PROBABILITY DENSITY FUNCTION ... 30 

C-3 Y-AXIS IMU ALIGNMENT ERROR PROBABILITY DENSITY FUNCTION ... 31 

C-4 Z-AXIS IMU ALIGNMENT ERROR PROBABILITY DENSITY FUNCTION ... 32 


v 



1.0 SUMMARY 


A comprehensive accuracy analysis of orbit IMU alignments using the Shuttle 
star trackers has been completed and the results are presented herein. Monte 
Carlo techniques were used in a computer simulation of the IMU alignment hard- 
ware and software systems to: (1) determine the expected STS-1 (manual mode) 

IMU alignment accuracy, (2) investigate the accuracy of alignments in later 
Shuttle flights when the automatic mode of star acquisition may be used, and 
(3) verify that an analytical model previously used for estimating the align- 
ment error is a valid model. 

In summary, the analysis results do not differ significantly from expecta- 
tions. The standard deviation in the IMU alignment error for STS-1 alignments 
was determined to be 68 arc seconds per axis. This corresponds to a 99.7% 
probability that the magnitude of the total alignment error is less than 
258 arc seconds. 



2.0 INTRODUCTION 

IMU alignments are performed in order to reposition each IMU inertial plat- 
form to a desired orientation with respect to the Mean of 1950 coordinate sys 
tern. The platforms must be realigned periodically because they do not remain 
perfectly inertial but drift away from tneir desired orientations. Figure 
2.0-1 illustrates the various coordinate systems associated with the IMU 
alignment. Each circle represents a coordinate system and the connecting 
lines represent coordinate transformations. In order to reposition a plat- 
form, its present orientation must first be determined. This can be accom- 
plished by measuring the positions of two stars relative . to the platform. The 
star measurements are nominally acquired by using the star trackers (ST). The 
Mean-of-1950 to measured IMU platform transformation matrix, M, is deter- 
mined from these two star measurements (Reference 1). The measured 
transformation matrix, M, is different from the actual transformation 
matrix. A, due to the star tracker and IMU measurement errors. Next, 
a platform repositioning matrix, T, is computed from 

T * RM t (1) 

where the matrix R defines the transformation from the Mean-of-1950 coordi- 
nate system to the desired IMU platform orientation. Finally, IMU torquing 
commands are extracted from this matrix and applied to the IMU. Note that 
the torquing angles are applied to the actual IMU; therefore, toe final orion 
tation of the IMU will not coincide with the desired orientation. The angu- 
lar displacement between the repositioned IMU platform and the desired orien- 
tation is referred to as the IMU alignment error. Assuming there are no 
errors associated with the torquing of the IMU, the alignment error transfor- 
mation matrix, denoted by TET T in Figure 2.0-1, is similar to the transforma- 
tion matrix E. In other words, both transformations represent the same 
eigenaxis rotation in inertial space. The alignment error matrix E can be 
computed from 

E « MA t (2) 

assuming that the actual platform orientation matrix A is known. 

In this analysis, the IMUCAl simulation program was used to determine the 
properties of the IMU alignment error matrix E using Monte Carlo techniques. 
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3.0 DISCUSSION 

The IMUCAl simulation program has been modified to compile statistical data 
on I MU alignment accuracy for multiple sample cases. For each alignment sam- 
ple (one Monte Carlo cycle), an alignment error matrix E is computed from 
equation 2. The Mean-of-1950 to measured platform transformation matrix M is 
computed as an intermediate step in the onboard alignment process while the 
actual platform orientation matrix A is available from the environment. 

Once the error matrix E is determined, an alignment error vector * is then 
computed as follows: 


* * 2 — sin'l iQl 

IQI 


(3) 


where 

Q 


1 

2 Vl + Eli + E 22 + E 33 


123 ‘ 132 
E 31 - E 13 

Ul2 - 4l\ 


( 4 ) 


denotes the vector part of a quaternion equivalent to the alignment error 
matrix E. The vector ♦ defines the magnitude ard direction of the angular 
displacement of the measured IMU platform with respect to the actual platform 
orientation. Note that * is expressed in the actual platform coordinate sys- 
tem. For each case, the simulation program computes the mean and the stan- 
dard deviation of each of the alignment error vector components 4 > x » 4 >y, as 
well as the alignment error vector magnitude |<&|. These statistics were 
computed using the following equations (k = number of samples). 


MEAN (per axis, j * x, y, z) 

1 k 

U j * - 2 (4>j)i 

J k 1*1 J 

STANDARD DEVIATION (per axis, j r x, y, z) 

V I f k " I k 7* 

w U Wj ^ ' i <*1 Wj)l) . 

M EAN (total) 

1 k 


» * r * d*Di 

k i*l 


(5) 


( 6 ) 


(7) 


STANDARD DEVIATION (total) 


V n. r k , 1 k 

— [ i J l ii, (1*1)1 
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1) 2 ' 

- 


( 8 ) 
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A total of eighteen different IMU alignment cases were simulated. In each 
case, j.OO samples of the IMU alignment error * were generated for each- IMU. 
These cases, which are presented in Table 4.1-1, simulated various 
combinations of star pair spatial separations (<S) and temporal separations 
{ t ) . The STS-1 trajectory (cycle 2) was simulated in every case and the 
simultaneous right angle pair case (case 1) conformed to the first STS-1 
orbit alignment in detail. The selected alignment stars were Achernar and 
Alpheratz which were sighted . imUltaneously at 2:45 GET by the -Z and the -Y 
star trackers, respectively, 'he remaining cases simulated alignments for 
later flights when the automatic mode of star acquisition may be used. In 
the automatic mode, stars, of opportunity are acquired under software direction 
and a wide range of spatial and temporal separations are possible. Software 
alignment star selection limits the separation ranges to 

35° < 6 < 145° 

0 < t < 90 minutes 

The temporal separation limit (sighting age limit) is applied only when data 
for four stars are available. 

For this analysis, angular pair separations between 30 degrees and 90 degrees 
and temporal separations less than 90 minutes were investigated. Angular 
separations greater than 90 degrees were not simulated because the alignment 
error function is symmetric with respect to right angle separations (see 
Section 4.0). 

For each sample, the IMU, ST, and Navigation 8ase bias errors were 
initialized at random, based on current estimates of the hardware system 
performance (Appendix A). Initialization also included the positioning of 
the three IMU's in a random inertial orientation with a fixed relative 
skew. Each IMU had 2 hours and 45 minutes of accumulated bias drift error at 
the time of the first star sighting. In those cases where a finite temporal 
separation was simulated, the first star was sighted by the -Z star tracker 
at 2:45 GET. After the desired amount of time had elapsed (30-90 minutes), 
an attitude maneuver was executed to place the second star in the -Y star 
tracker field of view (FOV). In each case, the alignment was performed 
immediately after the sighting of the second star. The second star 
measurement, therefore, was not corrupted by IMU drift error. 
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4.0 RESULTS 


Oetailed simulation results for each analysis case are givan in Tables 
8-1 through B-6 in Appendix B. Since several alignment error indicators 
are used an explanation of each is given below and summarized in Table 
4.0-1. The mean and standard deviation are defined by equations 5-8. 

TABLE 4.0-1 IMU Alignment Accuracy Indicators 



PER AXIS 
INDICATORS 

3 AXIS OR TOTAL 
ERROR INDICATORS 

MEAN 


V 

STANDARD DEVIATION 


a 

RMS 

; H9|| 

0) 


4.1 THE RMS IMU ALIGNMENT ERROR 

The root mean square (RMS) indicator 

u = Ve(|*| 2 ) = Vu 2 To 2 (9) 

can be expressed in terms of the spatial and temporal star pair separations, 
the variance of the single axis star sighting measurement error <? 0 2 , and the 
variance of the single axis IMU platform drift rate o j 2 
as follows (see Appendix 0 for derivation). 

u) * "yfaj- (1 + 2 csc 2 $) + Oj 2 t 2 csc 2 ^ (10) 

As can be seen from this equation, an optimum value of the alignment error 
occurs for a spatial separation of 90 degrees and a temporal separation of 
zero. The RMS indicator Is suitable for estimating the IMU alignment accu- 
racy since in a single number, it includes a measure of both the mean error 
and the variance of the error. It is also valid for use in the comparative 
analysis of different cases of spatial and temporal star pair separations. 

An estimate of the standard deviation of the single axis star position meas- 
urement error, cu, for star tracker measurements (Reference 1) and the stan- 
dard deviation or the single axis IMU bias drift rate, o^, (Appendix A) are 
given by 

o Q * 71.5 sec 

■ 0.02 sec/sec 

These performance estimates In conjunction with equation 10 were used to 
compute the analytical estimates of the RMS IMU alignment error for the 18 
cases of spatial and temporal st?r pair separations listed in Table 4.1-1. In 
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TABLE 4.1-1 RMS TOTAL IMU ALIGNMENT ERROR 



RMS ALIGNMENT ERROR (sic) 



Case 

6 

■ 

Analytical 

Estimate 

Monte Carlo 
Result 

1 

90 

0 

124 

121 

2 

90 

30 

129 

130 

3 

90 

60 

143 

148 

4 

90 

90 

164 

163 

5 

75 

15 

128 

121 

6 

75 

45 

139 

141 

7 

75 

75 

157 

159 

8 

60 

0 

137 

131 

9 

60 

30 

143 

131 

10 

60 

60 

160 

165 

11 

60 

SO 

185 

184 

12 

45 

15 

162 

147 

13 

45 

45 

177 

171 

14 

45 

75 

204 

193 

15 

30 

0 

215 

188 

16 

30 

30 

226 

193 

17 

30 

60 

258 

250 

18 

30 

90 

304 

300 




1x3 

124 

139 

160 

123 

134 

153 

131 

137 

155 

181 

155 

171 

199 

205 

217 

251 

298 









order to validate this model, also included in the table are the IMUCAL simu- 
lation derived values of the RMS alignmer* error. Note that the analytical 
estimates compare very well with the Monte Carlo simulation results. This 
verifies that the model (equation 10), used for estimating the alignment 
error is a valid model. There is, therefore, a high degree of confidence in 
this model based on both this verification and the original, rigorous develop- 
ment of the model. Conversely, many simplifying assumptions were made to ap- 
proximate the star sighting error, o 0 ; therefore, there is still some uncer- 
tainty in the estimated value. In order to determine a more accurate esti- 
mate of cr 0 , a curve fit to the simulation results was made using equation 10 

and the known standard deviation of the simulated IMU drift rate (0.02 

sec/sec). The tesults of the curve fit yields the following value for 

a c = 68.4 sec 

This value compares favorably with the estimated value. Table 4.1-1 contains 
data for the best fit surface to the simulation results using this 
empirically derived value for o 0 . The RMS IMU alignment error (equation 10) 
is plotted in Figure 4.1-1 as a function of the spatial and temporal star 
pair separations using the simulation derived value of o 0 . On this plot, 
there are several time scales corresponding to different values of the 
standard deviation in the single axis IMU drift rate, o d . 


4.2 THE SINGLE AXIS IMU ALIGNMENT ERROR 

For the purpose of navigation analysis, it is often desired to know the align- 
ment accuracy with respect to the individual axes of an inertial coordinate 
system. The single-axis alignment error properties, defined by equations 0- 
4 and D-5, in Appendix 0, apply only when the error vector * is resolved into 
the star pair relative coordinate system. The star pair coordinate system, 
however, can assume any orientation with respect to either the Mean-of-1950 
coordinate system or the IMU platform coordinate system. If the orientation 
of the star pair coordinate system is assumed to be uniformly random in 
three-space, the alignment error will be isotropic In an inertial ref- 
erence coordinate system. The distributions of the components, however, 
will not in general be Gaussian and, furthermore, cannot be solved for 
in closed *orm. Two important properties of this distribution can be 
Inferred from the symmetry of the star pair coordinate system components 
and the uniformity of the coordinate system's random orientation: 

U x * Uy * U 2 ■ 0, and (l?) 

a x ■ Oy » o 2 » ui 0 . (12) 

The subscripts x, y, and z are redefined to denote the coordinate axes of the 
inertial reference coordinate system. The equation for the RMS IMU alignment 
error (equation D-12) is valid for any coordinate system; therefore, equations 
D-12 and 12 can be combined to yield a solution for the variance of the 
single-axis IMU alignment error. 

a 0)2/3 (13) 
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FIGURE 4.1-1 RMS IMU ALIGNMENT ERROR 
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The standard deviation of the single axis I MU alignment error is plotted 
as a function of the spatial and temporal star pair separations in Figure 
4.2-1. The empirically derived value for o Q was used to generate the 
plot. There are several time scales corresponding to different values 
of the standard deviation in the single-axis IMU drift rate. 


The assumption concerning the isotropic nature of the IMU alignment error in 
the IMU platform coordinate system requires closer examination. The total 
alignment error is a vector sum of all the system measurement error components. 
Each error component has an associated direction that is determined by 
the relative orientation of the star tracker sensor measurement plane 
with respect to the IMU platform. Since the IMU platform was initialized 
at a random position in each Monte Carlo cycle, the distributions of the 
error components are expected to be isotropic. To test this assumption, 
an error covariance matrix was computed for the STS- 1 case (manual ST 
mode), using all 300 samples (that is, data from the 3 IMU's for Case 
1 ): 


4369 

161 

229 

161 

4416 

60 

229 

60 

5251 


(14) 


This covariance matrix is expressed in the IMU platform coordinate system and 
the units of the elements are sec . The covariance matrix defines an error 
ellipsoid by the equation 


K 2 = [x Y z] 



(15) 


The probability that an alignment error sample ^ will be inside the 
ellipsoid is a function of the parameter K. The lengths of the principal 
axes of the ellipsoid are given by 

ojj K * 73 K sec 

oy K * 67 K sec 

o| K » 65 K sec 

Since the probability ellipsoid approximates a sphere, the assumption of 
equal single-axis variances is reasonable. 



(16) 


The equation of the probability ellipsoid (equation 15), therefore, 
simplifies to that of a sphere with radius cj 0 K. 

(<d 0 K) 2 « X 2 + Y 2 .+ Z 2 (17) 
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FIGURE 4.2-1 

STANDARD DEVIATION OF THE SINGLE AXIS I MU ALIGNMENT ERROR 
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Note that the semiaxis (o^K), which is approximately parallel to the Z 
platform axis, is slightly longer than the other semiaxes. Covariance 
matrix diagonal elements were computed for the remaining 17 cases (using 
all 300 samples as in case 1) and, in all but one case (case 18), the 
variance in the Z-axis error was slightly larger than the variance in 
both the X- and Y-axis error. The variance in the X-and Y-axis errors 
were approximately equal in each case. This characteristic is a result 
of the choice of coordinate system. The azimuth resolver measurement 
error is not isotropic, but is unique since it is always parallel to the 
Z platform axis. The alignment error ellipsoid, therefore, is slightly 
elongated in the direction of the azimuth platform axis. For navigation 
simulations using the comprehensive IMU hardware model described in the 
Onboard Navigation Systems Characteristics document, the following values 
of the one sigma IMU alignment errors about the platform axes for $TS- 
1 alignments (6 * 90° 4 t = 0) are recommended: 

... 

o x * b6 sec 
oy * 66 sec 
o z * 72 sec 


4.3 THE MAXIMUM EXPECTED IMU ALIGNMENT ERROR 

It is often useful to know the range of the alignment error magnitude. Worst 
case is usually assumed to be three standard deviations from the mean for sin- 
gle dimensional random variables. There ’s a 99.74% probability that the sin- 
gle dimensional normally distributed random variable W is within three stan- 
dard deviations of the mean. 

P( u -3o < W < y + 3o) * 0.9974 

In the three dimensional case, the range associated with a given probability 
is an ellipsoidal volume centered about the mean. The mean IMU alignment 
error is zero for all three platform axes (equation 11) and the range 
is approximated by a spherical volume of radius u>qK (equation 17). The 
probability that an alignment error sample is contained within this volume 
is a function of the distribution and the value of K. A worst case IMU 
alignment error sample is defined to have a length that is equal to the 
radius of a sphere t.iat encompasses 99.74% of the population of alignment* 
error samples. If the alignment error ♦ is assumed normally distributed , 
then the value of K for this sphere is 3.77. 


* The IMU alignment error is normally distributed only in the STS-1 (manual 
ST mode) case (6 a 90° 4 t a 0 min). It has been shown, however, that the 
actual distribution for the general case closely approximates the normal 
distribution for the entire range of star pair temporal separations and 
a limited range of spatial separations. 


60° < 6 < 120° 

0 min < t < 90 min 
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F(j*| < 3.77 u 0 ) = 0.9974 

(18) 

Therefore, 


l*lmax * 3 * 77 “o 

(19) 


Note that the worst case error is a total error and that the single axis com- 
ponents will be smaller. Lines of constant maximum IMU alignment error 
(equation 19) are plotted as a function of the star pair temporal and spatial 
separations in Figure 4.3-1. In spite o* the assumptions made, the estimated 
maximum IMU alignment error as defined by equation 19 compares favorably with 
the simulation results in Appendix 8. 


4.4 THE STS-1 CASE 


The following sunmarizes the statistics for the STS-1 type (manual ST mode) 
alignments (see Table 4.0-1 for nomenclature). Derivation of these statistics 
is presented in Appendix 0. 

U x * Uy * Uz ~ 0 

o x * Oy 3 o z = 68 sec 
u = 109 sec 
o * 46 sec 


iD 0 3 68 sec 
i*Lax 3 258 


4.5 SPECIAL CASES 

The Eighteen test cases in Table 4.1-1 simulated star pair separations that 
ware less than or equal to 90 degrees. These cases were limited to this 
range because cf the symmetry of the alignment error function (equation 
10)'with respect to 90 degrees. An additional case, however, was simulated 
to verify the property of symmetry. Table 4.5-1 summarizes the results 
of this case. 
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FIGURE 4.3-1 

ESTIMATED MAXIMUM TOTAL IMU ALIGNMENT ERROR 
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TABLE 4.5-1 RMS ALIGNMENT ERROR - SPECIAL CASES 



RMS ALIGNMENT ERROR (sTc) 



Analytical 

Monte Carlo 

Curve 

Case 

6 t 

Estimate 

Result 

Fit 

o> 

V 

CO 

O 

O 

150 60 

258 

241 

251 

post-sleep 

90 0 

124 

125 

118 


One final case was simulated to assess the feasibility of star acquisition 
with large amounts of IMU drifc. The longest period between alignments 
for STS-1 will occur during the crew sleep periods. The first post-sleep 
alignment for STS-1 was simulated for E00 Monte Carlo cycles. This align- 
ment occurs ten hours and thirty minutes after the previous alignment. As 
in the other cases, the standard deviation in the per axis IMU drift rate 
was 0.02 degrees per hour. In the 200 attempts to acquire both stars in 
the offset mode, 28 failed. This corresponds to an 86% probability of 
success for offset mode acquisition when performing post-sleep alignments. 
This probability is based on the simulation of a perfect star tracker and 
a celestial sphere devoid of everything but navigation stars. Hardware 
anomalies and non-navigation trackable objects would decrease this estimate 
slightly for actual flight. The simulation was configured to investigate 
acquisition probability for only the offset mode. When the desired star 
is rot found in the offset node, the tracker is commanded to search the 
full field of view. The simulation program, however, cannot be used to 
determine the acquisition probability for full field search mode because 
of uncertainties in debris density and angular rates. 




5.0 CONCLUSIONS AND RECOMMENDATIONS 


In suimiary, the orbit IMU alignment accuracy analysis results compare favor- 
ably with previous estimates. The simulation-derived value for the standard 
deviation in the single-axis IMU alignment error (uu) for STS-I alignments is 
68 arc seconds. This corresponds to a 99. 7% probability that the magnitude 
of the total alignment error is less than 258 arc seconds. The. simulation 
program also* verified the analytic algorithm for estimating the IMU alignment 
accuracy (equation 10) as a function of the spatial and temporal separation 
of measured star pairs. 


5.1 CREW PROCEDURES 

Based upon this analysis, several recommendations concerning crew procedures 
for IMU alignments are made. Before the IMU platform is repositioned, the 
crew will check the displayed IMU torquing angles for reasonableness. If the 
angles are less than an acceptable limit, then the repositioning of the IMU 
platforms is executed, otherwise, the star measurements must be repeated. The 
torquing angle limit must be sized to include the initial alignment error, 
the IMJ drift error since the initial alignment, and the final alignment 
error. The longest intervals between alignments are during the crew sleep pe- 
riods and can be as long as 10.5 hours for STS-1. The variance in the 
torquing angles, therefore, is 

a T 2 * wo 2 + (10.5 ad) 2 + too 2 
a 762 sec 

The maximum expected magnitude of the torquing vector t is defined by 
P(|t| < 3.77 a T ) * 0.9974 
Mmx • 0-8° 

For simplicity, present crew procedures require that the largest torquing 
angle component be less than the acceptable limit; therefore, the recommended 
torquing angle limit criteria for alignment is 

t x < 0.8° 

x y < 0.8° 

x 2 < 0.8° 

This recommended value confirms the value currently in the STS-1 crew 
procedures. This simplified torquing angle component test is crude and 
not as exact as a magnitude (|x|) test. This could be resolved by displaying 
| x | on the IMU alignment display. The torquing vector magnitude Is presently 
available in the flight software as the variable C6MV QANG in the procedure 
GX1_MIS_ANG. 

The displayed torquing angles are also crew evaluated during IMU alignment 
verification. After each STS-1 alignment, a second alignment is performed 


16 


to verify the first. If any one of the displayed verification torquing 
angles exceeds an acceptable limit, then the alignment and verification 
process must be repeated. The verification torquing angle limit must 
be sized to include the initial IMU alignment error and the verification 
alignment error. The variance in the verification angles will be 

°t 2 s u o 2 + “o 2 
a T = 97 sec 

The maximum expected magnitude of the verification torquing vector t is 
defined by 

P ( | t | < 3.77 o T ) = 0.9974 
Mmax = 

As in the alignment case, the crew performs a rough check by simply verifying 
that each torquing angle component is less than the torquing angle limit. 

The recommended verification criteria, therefore, is 

T„ < 0.1° 

T, < 0.1° 

Tj < 0.1° 

Again, this recommended value confirms the value currently In the STS-1 
craw procedures. 
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APPENDIX A 


STAR TRACKER, IMU, AND NAV BASE 
ERROR MODEL PARAMETERS 



Error Source 


IMU ERROR MODEL 


la per axis error 
(arc sec) 


Resolver errors 

Resolver bias 30 

Random noise 12 

Quantization 11.5 

Sinusoidal bias 

1st harmonic 7.6 

8th harmonic 19.0 

9th harmonic 4.2 

16th harmonic 20.0 

Gimbal nonorthogonalities 

Pitch-to-outer roll 30 

Outer roll to case 20 

Gyro errors 

Bias drift .02/sec 


NAV BASE ERROR MODEL 


Error Source 


IMU LRU installation 
Thermal nav base bending (ST to IMU) 
IMU to ST structural uncertainty 
ST LRU installation 


la per axis error 
(arc sec) 

20 

13.7 

7.1 

20 


STAR TRACKER ERROR MODEL 


Error Source lo per axis error 

(arc sec) 

Bias 42.4 

Random 10.6 
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APPENOIX B 

MONTE CARLO SIMULATION DERIVED 
IMU ALIGNMENT ERROR STATISTICS 
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TABLE B-l KEY TO SIMULATION PROGRAM OUTPUTS 



2; 



TABLE 8-2 IMU ALIGNMENT ERRUR STATISTICS FOR 90° STAR PAIR ANGULAR SEPARATIONS 
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TABLE B-3 IMU ALIGNMENT ERROR STATISTICS FOR 75° STAR PAIR ANGULAR SEPARATIONS 
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TABLE B-4 IMU ALIGNMENT ERROR STATISTICS FOR 60 u STAR PAIR ANGULAR SEPARATIONS 
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TABLE B-5 IMU ALIGNMENT EI^ROR STATISTICS FOR 45 STAR PAIR ANGULAR SEPARATIONS 
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0.0 THE RMS IMU ALIGNMENT ERROR 


The RMS IMU alignment, error indicator is an analytical tool that can be used 
for estimation of the IMU alignment accuracy. The RMS indicator was first 
introduced in Reference 2 and it is the purpose of this appendix to present 
the derivation of this alignment error indicator. 

0.1 DERIVATION OF ALIGNMENT ERROR INDICATOR 

Recall that the alignment error is completely defined by the rotation vector 
*. It was shown in Reference 2 that with the proper choice of coordinate 
system, the components of $ can be determined as a function of the alignment 
star pair angular separation 6, the age of the most recently sighted star t s , 
the age of the oldest star sighting t t , the star direction measurement errors 
S±» S||* T ll> and the IMU drift rates d x , dy, d z . 

4> x = " T 1 + (ts + tt) d x sin *s5 + (t s - tf) dy cos %5)/sin *s6 

4>y = ^(Sx + Tjl + (t s - t t ) d x sin + (t s + t t ) d y cos H<S)/cos HS (0-1) 

4>z 3 $11 + ts d z 

The coordinate system used is fixed relative to the alignment star direction 
vectors S and T. These vectors define a plane referred to as the pair plane. 
The X-axis of this coordinate system lies in the pair plane and bisects the 
acute angle between S and T. The Z-axis is perpendicular to the pair plane 
(positive toward SxT) and the Y-axis completes the right angle triad. Sx and 
Sm are the star measurement angular error components perpendicular and par- 
allel to the pair plane, respectively, for the most recently sighted star S. 
Both Sx and S|| are normal to the star direction S. Similarly, and Tj| 
are the star measurement angular error components for the oldest star T. 
Finally, d x , dy, and d z are the components of the IMU drift rate about the 
coordinate axes X, Y, and Z, respectively. 

In the automatic star acquisition mode, MJ alignments will usually be per- 
formed immediately following the acquisition of the second star. Hence, this 
is the case that is of most interest. If the age of the most recently 

sighted star is negligible (t s * 0 and t t * t), then equation 0-1 simplifies 

to: 

<t>x * ^(Sj. * !l + * d x sin - t dy cos *s6)/sin »*S 

$y s *j(Sx + lx - t d x sin H6 + t dy cos Hfij/cos *s6 (0-2) 

4>z 3 s ll 

Tie parameters 5 and t are referred to as the spatial and temporal separa- 
tions of the alignment stars, respectively and are predetermined for a given 
pair of star measurements. The system errors $x, $||, Tx, Tj|, d x , d y , and d 2 , 
on t'.ie other hand, are random variables which cause the alignment error vec- 
tor components to be random variables also. The star direction measurement 
errors, $x» Sjj, T^, and Tjj, are all assumed to be normally distributed, inde- 
pendent random variables with zero means. Since these errors are the sum of 
a large number of randan components (Reference 1), the Central Limit Theorm 
makes this a reasonable assumption. The IMU drift rates, d x , dy, d z , are 
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also assumed to be normally distributed, independent random variables 
with zero means. Additionally, the measurement errors and the drift rates 
are assumed to be isotropic in three space. The expectations and the 
variances of these random variables are summarized by the following equations: 


E(S a ) = E(S„) * E(T X ) * E(T„) - 0 

V(S A ) = V(S„) - VtT A ) * V(T„) * cr 0 2 

E(d x ) • E(dy) * E(d z ) = 0 

V(d x ) * V(d y ) * V(d z ) = a d 2 


(D-3) 


Given the aforementioned assumptions and the properties defined by equations 
D-2 and D-3, it can be shown that the alignment error components are 
also independent and normally distributed with expectations and variances 
defined as follows: 

E( 4>x ) * Ux * 0 

EUy) * Wy 3 0 (D-4) 

E(d>z ) a v»z * 0 

V(<D X ) s a x 2 s H(a 0 2 + H t 2 a d 2 )/sin 2 

V(<j>y) = ay 2 = H(o 0 2 + H t 2 a d 2 )/cos 2 (D-5) 

V(d> 2 ) = a z 2 = o 0 2 


The joint probability density function (pdf) for the alignment error components, 
therefore, is defined by: - 

P 2 M.Z M.l\ 


f(*x» $y» ♦ z) 


- H 


(2w)3/2 


a x°y°z 


4>x f - <J>y‘ <t>z‘ 


»x‘ 


°z 


( 0 * 6 ) 


Although this function completely defines the properties of the alignment 
error vector *, it is not very useful for several reasons. First of all, 
this function describes the alignment error properties in a coordinate 
system fixed relative to the star directions. The alignment error properties 
are most useful if they can be related to the IMU platform coordinate 
system or the Mean-of-1950 Shuttle reference inertial coordinate system. 

These two systems can assume any orientation relative to the alignment 
star pair. The directional characteristics of the alignment error vector, 
however, are of secondary importance. Of primary interest is the magnitude 
of the alignment error vector whi ch is given by 

i*| * V^x 2 + 4>y 2 + *z 2 (0-7) 

The pdf g(|$|) is the panacea for alignment error estimation because it 
is independent of the choice of coordinate system. Once the pdf g(|*|) 
is known, then the expectation and the variance of |4| can be determined. 
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Unfortunately, there is no closed form solution .or g(|$|). 

The RMS IM'J alignment error indicator is a serviceable alternative to 
this dilemma. Recall that the expectation and the variance of a random 
variable are related by 

E(|*| 2 ) » V(j*|) + E( |*|)2 (0-8) 

In the adopted notation, this is equivalent to 

, g2 + jj2 ( D-9 ) 

Although u and o cannot be solved for individually, ie sum of u 2 and 
cr can be solved for. The expectations of the sauares of both sides of 
equation D-7 is 

E(|*| 2 ) = E(<J> X 2 + <P y 2 + h 2 ) 

E(l*l 2 ) = E(<D X 2 ) + E(* y 2 ) + E(4> z 2 ) (0-10) 

E(|<t| 2 ) » V(D X ) + E($ x ) 2 + V(4* y ) + E($ y ) 2 + V(<D Z ) + E(d> z ) 2 

Since the error components have zero means (equation D-4), equation U-10 
simplifies to 

E( |$| 2 ) = V(4> x ) + V(<p y ) + V(<p z ) (D-ll) 

which in terms of the adopted notation can be written as 

uj 2 = a x 2 + Oy 2 + o y 2 (D-12) 

By combining equations D-5 and D-12, the RMS IMU alignment erro: , <u, 

can be expressed in terms of the soatial and temporal star pair separations 

as follows 

a) = -yo 0 2 (1 + 2 esc 2 6) + o d 2 t 2 esc 2 6 (D-13) 

D.2 THE STS-1 CASE 

Recall that in Section D.l it was stated that the pdf of the IMU alignment 
error magnitude g(|$|) could not be solved for in closed form. For a 
single specific case, g(|4>|) can be determined. When IMU alignments are 
performed during STS-1, the spatial and temporal star pair separations 
will be optimal. 

5 = 90° 

t = C min 

In this case, the variances of the alignment .-ror components (equation 
D-5) are simply 
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The component errors are normally distributed with equal variances and 
zero means. The pdf of |4>|, therefore, is defined by Maxwell's distribution. 



This function is plotted and superimposed on the simulation derived pdf 
in Appendix C. The probability density functions for the single axis 
alignment errors about the X, Y, and Z platform axes are also included 
in Appendix C. The expectation and the variance of the Maxwell distribution 
are defined by 


E(|*|) = 2a 0 V^ 

( D- 15 ) 
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